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ABSTRACT 

We present a catalog of 25 definite and 11 probable strong galaxy-galaxy gravitational lens sys- 
tems with lens redshifts 0.4 < z < 0.7, discovered spectroscopically by the presence of higher redshift 
emission-lines within the Baryon Oscillation Spectroscopic Survey (BOSS) of luminous galaxies, and 
confirmed with high-resolution Hubble Space Telescope {HST) images of 44 candidates. Our survey 
extends the methodology of the Sloan Lens ACS Survey (SLACS) to higher redshift. We describe 
the details of the BOSS spectroscopic candidate detections, our HST Advanced Camera for Surveys 
(ACS) image processing and analysis methods, and our strong gravitational lens modeling procedure. 
We report BOSS spectroscopic parameters and ACS photometric parameters for all candidates, and 
mass-distribution parameters for the best-fit singular isothermal ellipsoid models of definite lenses. 
Our sample to date was selected using only the first six months of BOSS survey- quality spectro- 
scopic data. The full five-year BOSS database should produce a sample of several hundred strong 
galaxy-galaxy lenses and in combination with SLACS lenses at lower redshift, strongly constrain the 
redshift evolution of the structure of elliptical, bulge-dominated galaxies as a function of luminosity, 
stellar mass, and rest-frame color, thereby providing a powerful test for competing theories of galaxy 
formation and evolution. 

Keywords: galaxies: distances and redshifts - galaxies: evolution - galaxies: high-redshift - gravita- 
tional lensing: strong - cosmology: large-scale structure of universe - surveys 



1. INTRODUCTION 

Strong galaxy-galaxy lensing produces unmistakably 
distorted, amplified, and multiplied images of a distant 
galaxy through the gravitational field of a foreground- 
galaxy. This phenomenon provides a direct astrometric 
measurement of the total mass distribution within the 
interior of the lens galaxy, including luminous and dark 
components. Strong lensing on galactic scales was first 
suggested by Zwicky (1937) and first discovered by Walsh 
et al. (1979). Since then, a number of dedicated and 
serendipitous surveys have produced a combined sample 
of over 200 galaxy-scale lens systems (for a recent re- 
view, see Treu 2010). Of the various survey methods, the 
most prolific means of identifying strong galaxy-galaxy 
lens candidates is spectroscopic discovery, which provides 
evidence for two galaxies along one line of sight, along 
with precise measurements of the lens and source red- 
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* Based on observations made with the NASA/ESA Hubble 
Space Telescope^ obtained at the Space Telescope Science Insti- 
tute, which is operated by the Association of Universities for Re- 
search in Astronomy, Inc., under NASA contract NAS 5-26555. 
These observations are associated with program 12209. Based 
on spectroscopic data from the Baryon Oscillation Spectroscopic 
Survey of the Sloan Digital Sky Survey HI. 



shifts (Warren et al. 1996; Bolton et al. 2004). Images 
taken with the Hubble Space Telescope {HST) or other 
high-resolution facilities then provide a confirmation of 
the lensing hypothesis and make precise (milliarcsecond) 
measurements of the relative positions of the compo- 
nents. This method has been employed by the Sloan Lens 
ACS (SLACS) survey (Bolton et al. 2006, 2008) to pro- 
duce a sample of 85 grade-A confirmed strong lenses and 
an additional 13 grade-B probable lenses (Auger et al. 
2009). 

Strong lensing is the best probe of the distribution 
of mass in the transition region between the baryon- 
dominated central regions of galaxies and the dark mat- 
ter halo (e.g., Kochanek 2005). By combining strong- 
lensing aperture masses with lens-galaxy luminosities, 
sizes, and velocity dispersions, the SLACS collaboration 
has arrived at the following results (Treu et al. 2006; 
Koopmans et al. 2006; Gavazzi et al. 2007; Bolton et al. 
2008; Koopmans et al. 2009, see also Jiang & Kochanek 
2007). 

1. The average total mass-density profile of early- 
type galaxies is nearly isothermal, corresponding to 
nearly flat galaxy rotation curves within the Ein- 
stein radius. 
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2. Models for early- type galaxies in which light lin- 
early traces mass are falsified at > 99.9% confi- 
dence even within the optical half-light radius. 

3. The average weak-lensing signal out to a few hun- 
dred kiloparsecs around the SLAGS lenses is consis- 
tent with an isothermal extrapolation of the strong- 
lensing constraints from the inner few kiloparsecs. 

4. The average density structure of elliptical galaxies 
is universal and does not vary systematically with 
mass and exhibits a small but significant intrinsic 
scatter that is not significantly correlated with any 
other observable quantity. 

5. The relationship between the lensing mass and the 
dynamical mass of SLAGS lenses has a logarithmic 
slope of LOS ±0.04. 

6. The variation in the dark matter fraction suggests 
a bulge-halo conspiracy^ with the combined rota- 
tion curve of both components together being flat- 
ter than that of either component independently. 

Many of these results were also found in earlier but 
smaller studies of lensed quasars (see e.g., Treu & Koop- 
mans 2004, or Rusin & Kochanek 2005, and references 
therein) . 

In order to conduct the next generation of large- 
scale structure surveys, the Sloan Digital Sky Survey 
III (SDSS-III) collaboration upgraded the Sloan Digi- 
tal Sky Survey I (SDSS-I) optical spectrographs (York 
et al. 2000) and initiated the Baryon Oscillation Spec- 
troscopic Survey (BOSS; Eisenstein et al. 2011). The 
significant hardware changes relative to SDSS-I are (1) 
an increase in the number of fibers from 640 to 1000, 
(2) smaller fiber diameters that subtend 2^' rather than 
3^' on the sky, (3) volume-phase holographic gratings 
with increased efficiency, and (4) larger and more sen- 
sitive GGDs that increase throughput and wavelength 
coverage relative to SDSS-I. With this new configura- 
tion, BOSS is obtaining spectra of approximately 1.5 
million luminous galaxies out to redshift z 0.7 from 
mid-2009 to mid-2014. This sample complements the 
SDSS-I luminous red galaxy (LRG) dataset that extends 
out to redshift z 0.35. LRGs are luminous, ellip- 
tical, bulge-dominated galaxies, with remarkably uni- 
form spectral energy distributions that are character- 
ized by a strong break at 4000 A (Eisenstein et al. 2001, 
2003). The BOSS luminous galaxy sample is composed 
of a low-redshift (LOZ) sample designed to select LRGs 
with 0.2 ^ z ^ 0.4 and a constant mass (GMASS) 
sample designed to select luminous galaxies with stel- 
lar masses in the range lO^^M© < < lO^^M© with 
0-4 < ^ < 0.7 (Eisenstein et al. 2011). Masters et al. 
(2011) performed visual inspection of 240 BOSS targets 
in the GOSMic Origins Survey field and determined that 
^ 75% of BOSS galaxies are "early- type" , with the re- 
mainder "late-type", and ~ 20% of BOSS galaxies are 
multiple systems. The combined sample will enable a 
baryon acoustic feature measurement (Eisenstein & Hu 
1998; Eisenstein et al. 2005) of the distance-redshift re- 
lation in several redshift bins. 

As a consequence of this design and survey strategy, 
BOSS presents an opportunity to discover a sample of 
strong galaxy-galaxy lenses at significant cosmological 
look-back time that is of comparable size and homogene- 
ity to the SLAGS sample at lower redshift. The two 



samples in combination may allow a unique study of the 
evolution of the structure of mass density and its de- 
pendence on other parameters such as stellar mass and 
stellar velocity dispersion. 

A complementary sample of higher-redshift lenses dis- 
covered in the Ganada-France-Hawaii Telescope Strong 
Lens Legacy Survey (SL2S; Gabanac et al. 2007; Ruff 
et al. 2011) suggests evolution in the slope of the mean- 
density profile to shallower slopes at higher redshifts. 
The key complementarity between SL2S and BOSS 
strong-lens surveys is that the former is more uniformly 
sensitive to selecting lenses across a range of angular 
scales (and hence mass scales), whereas the latter guar- 
antees the spectroscopic redshift completeness for both 
lensing and lensed galaxies that is necessary for an accu- 
rate conversion from observed angular scales to physical 
mass scales. 

This paper (Paper I) is the first in a series presenting 
the initial catalog of 45 spectroscopically selected strong 
galaxy lens candidates discovered from the first partial 
year of BOSS data. We refer to this survey as the BOSS 
Emission-Line Lens Survey (BELLS). Of the initial 45 
candidates, 44 were successfully observed with the Wide 
Field Gamera (WFG) of the HST Advanced Gamera for 
Surveys under Program ID GO-12209 (PI: A. Bolton), re- 
sulting in 25 confirmed strong gravitational lenses, along 
with an additional 11 possible lenses with either very 
faint candidate lensed features or significant dust fea- 
tures complicating the analysis. The remaining target, 
SDSS J0212+0027, is scheduled for HST visit. We show 
that the BELLS and SLAGS lenses trace a single evolv- 
ing massive galaxy population in terms of stellar mass, 
thus confirming the suitability of the combined sample 
to constrain the structural and dynamical evolution of 
massive galaxies down to the present-day universe. 

This paper is organized as follows. The spectroscopic 
discovery method is described in §2. The HST imaging 
data, reduction method, and the sample of 45 candidates 
are presented in §3, and the photometric analysis includ- 
ing the B-spline fits used for lens-galaxy subtraction and 
the de Vaucouleurs (1948) models used for measuring 
the magnitudes and half-light radii are presented in §4. 
We present our strong-lens modeling methods and re- 
sults in §5, employing singular isothermal ellipsoid (SIE) 
mass models in combination with both parametric and 
pixelized-grid source-plane surface-brightness models. A 
summary of results is presented in §6. Throughout this 
paper, we assume a standard general-relativistic cosmol- 
ogy with parameters (/i, Qm^ ^a) = (0.7, 0.3, 0.7) (Larson 
et al. 2011). 

2. SPEGTROSGOPIG GANDIDATE SELEGTION 

The premise behind the spectroscopic strong-lens can- 
didate selection of Bolton et al. (2004) was to search 
for multiple background emission-lines within the 3^' di- 
ameter solid angle covered by the SDSS-I spectroscopic 
fiber in the residual spectra found after subtracting 
best-fit galaxy templates to the foreground-galaxy spec- 
trum (also see Warren et al. 1996; Hewett et al. 2000; 
Willis et al. 2005, 2006). For our current search within 
the BOSS database, we apply much the same method, 
which is also described in Bolton et al. (2006, 2008), to 
the 2^' diameter solid angle covered by the SDSS-III spec- 
troscopic fiber. As a result of the blurring effect of at- 
mospheric seeing, spectroscopic lens selection is efficient 
even at Einstein radii of Oe ^ 2" (Arneson et al. 2011). 
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TABLE 1 
Emission-Line Wavelengths 



Emission Restframe SDSS-I SDSS-III 

Line Wavelength (A) Zma^ Zma^ 

(1) (2) (3) (4) 



[Oii]AA3727 3729 88 ^'^^ ^'^^ 

m 4102.89 1.22 1.39 

H7 4341.68 1.10 1.26 

H/3 4862.68 0.87 1.02 

O III] 4959 4960.30 0.83 0.98 

III] 5007 5008.24 0.82 0.96 

II] 6548 6549.86 0.39 0.50 

Ha 6564.61 0.39 0.49 

[Nil] 6583 6585.27 0.38 0.49 

[S II] 6716 6718.29 0.35 0.46 

[S II] 6730 6732.68 0.35 0.46 



Note. — Column 1 lists the emission-lines used in 



the spectroscopic discovery of background emission- 
line sources and Column 2 shows the rest-frame vac- 
uum wavelength of the emission-line. Columns 3 and 4 
provide the maximum redshift of the emission-line de- 
tectable with the SDSS-I and SDSS-III spectrographs, 
respectively. 



The full list of possible background emission-lines for 
which we search is given in Table 1, along with the cor- 
responding maximum redshift to which they can be de- 
tected with either the SDSS-I or the more red-sensitive 
BOSS spectrographs. The spectra are based on SDSS-III 
imaging (Gunn et al. 2006, 1998), available in the eighth 
SDSS public Data Release (DR8) (SDSS-III collabora- 
tion: Aihara et al. 2011). 

We summarize our selection procedure here for com- 
pleteness, and to reflect minor changes in the implemen- 
tation relevant for BOSS. 

1. Rescale the spectrum noise- vector estimates based 
upon the residual fluctuations in sky-subtracted 
sky-fiber spectra on the plate. 

2. Select all objects both targeted and spectroscop- 
ically classified as galaxies, without any redshift 
warning flags. In the redshift analysis, only 
"galaxy" and "star" templates were allowed, since 
the inclusion of "qso" templates in the analysis 
of relatively low signal-to- noise ratio (S/R) BOSS 
galaxy spectra yields an unacceptable number of 
spurious classifications and redshift s. There were 
133,852 galaxy spectra scanned in total. 

3. Mask the positions of common emission-lines in the 
spectra of the target-galaxies, re-fit a model to the 
BOSS target galaxy continuum using a basis of 
seven principal component analysis (PGA) eigen- 
spectra generated as described in SDSS-III collabo- 
ration: Aihara et al. (2011) (versus four PGA eigen- 
spectra used in the redshift analysis) , and subtract 
this continuum model from the data. 

4. Scan the residual spectra automatically for at least 
two of the lines Ha, H/3, [Olll] 5007, [Olll] 4959, or 
[On] A A 3727 detected at 4cr significance or higher 
at a common redshift greater than that of the 
BOSS target. We refer to these as "multi-line" de- 
tections. There were 2898 multi-line hits. 

5. Scan the residual spectra automatically for a sin- 
gle line detected at 6cr significance or higher that 



would be at a higher redshift than the BOSS tar- 
get if identified as [Oll]AA3727. We refer to these 
as "single-line" detections. There were 8211 single- 
line hits. 

6. Bin the single-line detections by observed- frame 
wavelength, and prune those that are in highly 
overpopulated bins associated with night sky 
emission-line subtraction residuals. 

7. Bin the single-line detections by target-galaxy rest- 
frame wavelength, and prune those that are in 
highly overpopulated bins associated with imper- 
fect foreground-galaxy feature masking and sub- 
traction. There were 2137 single-line hits remain- 
ing after the cuts of Steps 6 and 7. 

8. Veto any single-line [Oll]AA3727 identification 
that is more plausibly explained as Ha, H/3, or 
[O III] 5007, based upon the appearance of sec- 
ondary lines at 3a significance or higher. There 
were 1624 single-line hits remaining. 

9. Remove any single-line detections also found in the 
multi-line search. There were 1268 single-line hits 
remaining. 

10. Remove any detections in which the background 
redshift is not greater than the foreground redshift 
by the margin, zi, — zs > 0.05. There were 1303 
multi-line hits and 741 single-line hits remaining. 

11. Visually inspect the reduced spectra of all remain- 
ing selected systems at the position of the candi- 
date background emission-lines, to remove obvious 
data-reduction/data-quality artifacts. There were 
1226 multi-line hits and 507 single-line hits remain- 
ing after this step. 

12. Inspect neighboring fiber spectra on each plate at 
the wavelength of the candidate background emis- 
sion to rule out cross talk from bright emission- 
lines in neighboring spectra, and spatially localized 
auroral emission not subtracted by the global sky 
spectrum model. 

13. Inspect the individual exposure spectra that were 
co-added to give the final detection spectrum, to 
ensure that the candidate background emission-line 
feature is not present in only a single frame (as 
would be expected for e.g., cosmic-ray hits). 

14. Inspect the raw spectroscopic GGD data for all ex- 
posures at the position of [O ll] AA 3727 emission for 
single-line candidates, to ensure that the detection 
is not related to cosmetic blemishes in the GGD. 

15. Gompute a total S/R for all detections as the 
quadrature sum of the individual S/R values (i.e., 
cr-values) of all lines in Table 1 detected at a sig- 
nificance of 3(7 or higher — including those which 
we do not employ as primary detection triggers in 
Step 4 — and rank candidates by this total S/R. 

In addition to ranking candidates by total S/R, we also 
select candidates to maximize lensing probability. Bolton 
et al. (2008) show empirically that the lensing probability 
of a candidate system is a strongly increasing function 
of the predicted Einstein radius of the system using a 
simple singular isothermal sphere model based on the 
foreground-galaxy velocity dispersion and the redshifts 
of both foreground- and background-galaxies. 
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TABLE 2 

BOSS Properties of BELLS Candidate Systems 



System Name 




Plate-MJD-Fiber 








zs 


i-band Magnitude 




ReS (") 


c^BOSS (kms 1) 


Sample 


(1) 




(2) 




(3) 




(4) 


(5) 




[b) 




(7) 


(8) 


bUbb JUi5iU7.o7+UU4yuy, 


,0 


ooUd-5518Z-Ud79 


0. 


,5171 


1 


.3636 


19.81 ± 0.05 





.89 + 0, 


.21 


01 _L QO 

219 ± 39 


CMASS 


dDdd JUz1z14.8U+UUz719, 


,1 


A^or" c /I TO o/^oo 

4z3o-55479-UoU3 


0. 


,5372 


1, 


.0122 


1 /I T 1 00 

19.47 ± 0.03 


1, 


.17 + 0, 


.19 


100 1 00 
189 ± 39 


CMASS 


dUdd JU747z4.iz+5U55o7, 


,5 


OCT"? cconc ncci 
oo77-55zU5-U551 


0. 


.4384 





.8983 


1 00 _L no 
18.92 ± 0.03 


1 


.24 + 0, 


.22 


QOO _L cir\ 

328 ± 60 


CMASS 


dUdd JU747o4. 75+444859, 


,3 


oo7d-5518d-U581 


0. 


,4366 





.8966 


1 /I _L n no 
18.84 ± 0.03 


2 


.87 + 0, 


.50 


281 ± 52 


CMASS 


dDdd JU75754.1U+4c51c55c5, 


,8 


3o7o-5518o-UUd3 


0. 


,5146 


1, 


.2165 


1010 1 /I 
19.13 ± 0.04 


2 


.51 + 0, 


.55 




CMASS 


dDdd JUoU1U5.c5U+47z749, 


,6 


o r" o A rza^Ar" 0001 

3o84-55z4o-U3Ul 


0. 


,4831 


1, 


.5181 


1 T 1 C 

19.87 ± 0.05 





.57 + 0, 


.17 


98 + 24 


CMASS 


bDbb JU8z13L).od+373c>c>U, 


,7 


oT/^o rro/^o o/^oo 

37d0-55zd8-0d99 


0. 


,5056 





.9705 


1 /I 1 00 
19.24 ± 0.02 





.93 + 0, 


.12 


000 1 00 
203 ± 23 


CMASS 


bDbb JU83L)49.7c>+511oc>1, 


,8 


3695-55212-0142 


0. 


,5301 


1 


.3317 


1 r 1 /I 

19.50 ± 0.04 


1 


.10 + 0, 


.22 


o/^ 1 

268 ± 36 


CMASS 


dDdd JU8c)7z7.UZ+49c57Uc5, 


,7 


o/^oT CTCTOoo oocrc; 

3o97-55z90-0z55 


0. 


,5513 


1, 


.1896 


1 c 1 000 

19.50 ± 0.03 


1, 


.24 + 0, 


.27 


000 1 A 

208 ± 34 


CMASS 


dUdd JU84U54.i9+5U5i5o, 


,1 


oo97-55z9U-U9Zl 


0. 


,5534 


1 


.4348 


1 n 1 _L n no 
19.81 ± 0.03 





.65 + 0, 


.17 


oo^^ _L QO 

206 ± 32 


CMASS 


dUdd JU84izi.74+5Ui7ZU, 


,4 


Q^oT ccoon nnoT 
oo97-55z9U-0U87 


0. 


,5542 


1 


.2881 


1 n CO _L OQ 

19.50 ± 0.03 





.79 + 0, 


.17 


176 + 25 


CMASS 


dDdd JU9151o.dU— UU55UU, 


,6 


37oo-55zl3-0o75 


0. 


,4022 


1, 


.1705 


t A C 1 00 

18.46 ± 0.02 


1, 


.98 + 0, 


.28 


221 + 16 


CMASS 


dDdd J(j941Uz.d9— U1U4UZ, 


,5 


oToo ccTO/i/i ocrio 

3782-55244-0513 


0. 


,4610 





.9115 


1 TO 1 /I 

19.79 ± 0.04 





.94 + 0, 


.25 


180 + 50 


CMASS 


dUdd JU944z7.47— Ui474z, 


,4 


QToo ccco/1/1 no^o 
o782-55244-02o8 


0. 


,5390 


1 


.1785 


1 TO _L O/l 

19.70 ± 0.04 


1 


.35 + 0, 


.46 


00 A _L Q A 

204 ± 34 


CMASS 


dDdd J1U1o58.z9— UzU8c5c5, 


,3 


OTO/l CTCTO/^O OOOT 

3784-55269-0237 


0. 


,4699 


1, 


.0034 


1 /I 1 /I 
19.34 ± 0.04 


1, 


.71 + 0, 


.47 


179 + 26 


CMASS 


bubb J 1113941. lb— U(Jl4z4, 


,5 


0000 ccrooo 0100 

3833-55290-0109 


0. 


,3849 





.9141 


1 TO 1 f\ r\ A 

19.72 ± 0.04 





.71 + 0, 


.17 


107 + 21 


CMASS 


OT^oo T1 1 1 TOT TO r> 1 oono 

bDbb J 111737. 7z— 01331)8, 


,9 


oToo r r /I 0000 

3788-55246-0930 


0. 


,4670 


1 


.2860 


1000 1 OT 

19.39 ± 0.07 


2 


.54 + 1, 


.24 


001 1 r 

221 ± 59 


CMASS 


bDbb J 115944. o3—L)L)U7z8, 


,2 


<~) A <~) rroTo 00/^0 

3843-55278-0069 


0. 


,5793t 


1 


.3457 


10/^0 1 000 
19.63 ± 0.03 





.99 + 0, 


.22 


165 + 41 


CMASS 


dDdd J1z15U4.44+UU47Zo, 


,0 


3846-55327-0559 


0. 


,6423 


1, 


.2970 


1 OT 1 00 

19.27 ± 0.03 


1, 


.42 + 0, 


.21 


262 + 45 


CMASS 


bDbb Jlzzll3.z8— UzzU37, 


,8 


i^i^i^i^ rroi ooTr 

3777-55210-0075 


0. 


,4061 


1 


.0108 


1 OT 1 00 

19.07 ± 0.02 





.64 + 0, 


.09 


/1 1 /I 
240 ± 24 


CMASS 


dUdd J 1zz151.9z+38Ud1U, 


,5 


3965-55302-0586 


0. 


,5348 


1 


.2844 


1 TT _L /I 

19.77 ± 0.04 





.93 + 0, 


.22 


1 T _L A 

187 ± 48 


CMASS 


O T^O O T100/10TOO r>0/l100 

dDdd J 1z34z7.99— Uz41z9, 


,6 


OTTO CrC01 OOTC 

3778-55213-0375 


0. 


,4900 


1, 


.0159 


19.18 ± 0.04 


1, 


.61 + 0, 


.46 


122 + 31 


CMASS 


bubb J 1318z9.39— UlU4zl, 


,6 


4004-55321-0121 


0. 


,6591 


1, 


.3959 


t r\ c c 1 /I 

19.66 ± 0.04 


1, 


.06 + 0, 


.31 


1 TT 1 OT 

177 ± 27 


CMASS 


OT^OO T1 ooTri 01 1 o/^or>i o 

bDbb J133751.31+3dZL)18, 


,1 


000/^ rrooo 0T1 

3986-55329-0713 


0. 


,5643 


1 


.1821 


1001 1 000 
19.01 ± 0.03 


1 


.60 + 0, 


.33 


225 + 35 


CMASS 


dDdd J1344z7.9o+3z58z4, 


,6 


ooc/^ ccro/^o 0000 

3856-55269-0309 


0. 


,4750 





.9777 


1 T 1 000 

19.57 ± 0.03 





.78 + 0, 


.13 


100 1 00 
198 ± 22 


CMASS 


o T^o o Tio/icrr>T/^o r>ioooo 

dDdd J1345L)7.o3— L)1z939, 


,8 


A r\ A A c o/^oo 

4044-55359-0603 


0. 


,4884 


1, 


.2992 


1000 1 /I 
19.92 ± 0.04 





.57 + 0, 


.13 


000 1 OCT 

200 ± 25 


CMASS 


CT^cc TiQ/inin on i q/^ioqo 
dUdd J 13491U.3U+3d1z39, 


,7 


QOCO CCO/IQ nooT 

3852-55243-0327 


0. 


,4396 





.8926 


1000 _L 00 

18.88 ± 0.02 


2 


.03 + 0, 


.31 


178 + 18 


CMASS 


dUdd J 135Zl8.99+3Zlo51, 


,8 


3861-55274-0393 


0. 


,4634 


1 


.0341 


1 00 _L OQ 

19.28 ± 0.03 


1 


.35 + 0, 


.33 


161 + 21 


CMASS 


OT^OO T1 /ICrOCO 1 /I 1 OOOO/IO 

dDdd J145z59.14+33z349, 


,4 


TO CTCTOTO OCTOCT 

3870-55273-0595 


0. 


,5602 


1, 


.1896 


1010 1 /I 
19.19 ± 0.04 


1, 


.53 + 0, 


.31 


OO/^ 1 A 

296 ± 42 


CMASS 


dUdd J 15U345.d9+3zz54z, 


,1 


QOT/^ cco/ic ncoc 
3876-55245-0585 


0. 


,6059 


1 


.4489 


10 A C\ _L r\Ci 

19.49 ± 0.06 


3 


.65 + 1, 


.18 


162 ± 68 


CMASS 


dUdd J 15zzU9.54+z91Uz1, 


,9 


QOTo ccco/1/1 nooo 
3879-55244-0299 


0. 


,5553 


1 


.3108 


1 T _L OQ 

19.67 ± 0.03 


1 


.08 + 0, 


.27 


1 _L OT 

166 ± 27 


CMASS 


OT^OO T1 CrOTOO OT 1 OOOOO/^ 

dDdd J15373U.z7+L)ZzU3o, 


,9 


Ar\nA c 0000 

4054-55358-0809 


0. 


,4825 


1, 


.0729 


10/^0 1 /I 
19.62 ± 0.04 





.57 + 0, 


.15 


199 + 25 


CMASS 


dDdd J154118.5o+181z35, 


,1 


OOOT C 00/^1 

3937-55352-0361 


0. 


,5603 


1, 


.1133 


1000 1 000 
19.82 ± 0.03 





.59 + 0, 


.00 


174 + 24 


CMASS 


SDSS J154246.33+162951, 


,8 


3932-55337-0331 


0. 


,3521 


1, 


.0233 


18.39 + 0.02 


1, 


.45 + 0, 


.13 


210 + 16 


LOZ 


SDSS J154503. 57+274805, 


,3 


3952-55330-0927 


0. 


,5218 


1, 


.2886 


18.90 + 0.03 


2 


.65 + 0, 


.47 


250 + 37 


CMASS 


SDSS J160113. 27+213833, 


,9 


3935-55326-0093 


0. 


,5435 


1, 


.4461 


19.75 + 0.03 





.63 + 0, 


.12 


207 + 36 


CMASS 


SDSS J161109.80+170526, 


,6 


4072-55362-0847 


0. 


,4766 


1, 


.2109 


19.66 + 0.04 


1, 


.33 + 


.38 


109 + 23 


CMASS 


SDSS J161523. 55+205636, 


,3 


4057-55357-0610 


0. 


,5894 


1, 


.3103 


18.42 + 0.03 


5 


.92 + 1, 


.03 


210 + 34 


CMASS 


SDSS J161846. 74+193026, 


,7 


4057-55357-0201 


0. 


,4906 


1, 


.0164 


19.93 + 0.04 





.76 + 0, 


.22 


275 + 53 


CMASS 


SDSS J162754.57+191004, 


,9 


4060-55359-0075 


0. 


,5251 


1, 


.4892 


19.59 + 0.03 





.59 + 0, 


.11 


236 + 33 


CMASS 


SDSS J163150.33+185404, 


,1 


4064-55366-0790 


0. 


,4081 


1, 


.0863 


17.70 + 0.01 


2 


.07 + 0, 


.10 


272 + 14 


LOZ 


SDSS J163714.58+143930, 


,1 


4065-55368-0515 


0. 


,3910 


0, 


.8744 


19.49 + 0.03 





.89 + 0, 


.15 


208 + 30 


CMASS 


SDSS J212252. 04+040935, 


,5 


4081-55365-0241 


0. 


,6261 


1, 


.4517 


19.29 + 0.03 


1, 


.76 + 0, 


.25 


324 + 56 


CMASS 


SDSS J212510.67+041131, 


,6 


4080-55471-0621 


0. 


,3632 


0, 


.9777 


18.21 + 0.01 


1, 


.47 + 0, 


.09 


247 + 17 


LOZ 


SDSS J230335. 17+003703, 


,2 


4207-55475-0933 


0. 


,4582 


0, 


.9363 


18.96 + 0.02 


1, 


.35 + 0, 


.23 


274 + 31 


CMASS 


Note. — Column 1 provides 


the SDSS System Name 


in 


terms of truncated J2000 R.A. and decl. in 


the format HHMMSS.ss+DDMMSS 


.s. Column 


2 provides the Plate-MJD-Fiber of the primary spectrum for the BOSS target (Eisenstein et al. 2011) 


a 


s of v5_4_45 of the BOSS pipeline. 


Columns 3 



and 4 provide the foreground BOSS spectroscopic redshift and the spectroscopically detected background emission-line source redshift, respectively. 
Columns 5 and 6 provide the BOSS de Vaucouleurs model i-band magnitudes and intermediate axis effective radii, respectively (Fukugita et al. 
1996; SDSS-III collaboration: Aihara et al. 2011). Column 7 provides the velocity dispersion computed by the BOSS pipeline (v5_4_45), uncorrected 
for aperture effects, using redshift-error marginalization and restricted stellar template sets as described in Shu et al. (2011). Column 8 indicates 
whether the candidate was selected from the BOSS LOZ or CMASS sample. All the BELLS candidates in this table were observed (or will be 
observed in the case of SDSS J0212+0027) by HST cycle 18 Program ID GO-12209. The HST-ACS /8i4-band properties are listed in Table 3, 
including our classification of morphology, multiplicity, and lens grade, with our justification of BELLS confirmed gravitational lenses in Table 4. 
t The lens redshift quoted for SDSS J1159-0007 differs from the BOSS v5_4_45 galaxy redshift that will be reported in the ninth SDSS pubhc Data 
Release (DR9) as z^noqso, due to the background emission-line flux. 

For BOSS lens candidates, we do not have precise ve- 
locity dispersion measurements to use in this calculation 
nor do we want to introduce a selection bias by using lu- 
minosity as a proxy. Hence, we assume a fiducial velocity 
dispersion of 250 kms~^ for all candidates and calculate 
a predicted "Einstein radius" given the foreground and 
background redshift s. This fiducial value is taken from 
the median of all well-measured velocity dispersions in 
the subset of the Bolton et al. (2008) SLAGS lens sam- 
ple having early-type morphology. We then compute the 
physical radius subtended by this predicted Einstein ra- 
dius at the redshift of the lens galaxy and require it to 
be greater than 4.5 kpc. Since no actual magnitude or 



velocity dispersion information is used to make this cut, 
it can only bias the sample in terms of foreground and 
background redshift distributions. The cut corresponds 
to a predicted Einstein radius of 0^.^84 at = 0.4, and 
0^63 at = 0.7. We verified that the sample selected 
by this cut does produce a reasonably uniform distribu- 
tion of candidates over the available redshift range in the 
parent sample. The fiducial values of = 250 kms~^ 
and 6e = 4.5 kpc are not especially significant and can 
be traded off against one another with no effect on the 
selected sample. Effectively, our cut serves to favor can- 
didates with higher background redshifts at a given fore- 
ground redshift, and pushes to slightly smaller angular 



BELLS I 



5 



Restframe Wavelength [A] 
3700 3725 3750 



Restframe Wavelength [A] 
3700 3725 3750 



Restframe Wavelength [A] 
3720 3740 4400 4600 4800 5000 





H7 



[Oiii] 5007"; 

[Oiii] 4959- 



8750 8800 8850 

Observed Wavelength [A] 

SDSS J0151+0049 

Restframe Wavelength [A] 
3720 3740 4840 4860 



7450 7500 7550 

Observed Wavelength [A] 

SDSS J0212+0027 

Restframe Wavelength [A] 
3700 3725 3750 



7100 8500 9000 

Observed Wavelength [A] 

SDSS J0747+5055 

Restframe Wavelength [A] 
3700 3725 3750 



- [On 


1 _ 

3727 - 


- H 

W-VWvA< 


13 ■ 


" , , 1 , , 










7050 7100 9200 9250 

Observed Wavelength [A] 

SDSS J0747+4448 

Restframe Wavelength [A] 
3700 3725 3750 3775 



8200 8250 8300 

Observed Wavelength [A] 

SDSS J0757+4313 

Restframe Wavelength [A] 
3700 3725 3750 



9350 9400 9450 

Observed Wavelength [A] 

SDSS J0801+4727 

Restframe Wavelength [A] 
3700 3725 3750 




" [On 


1 1 1 1 1 1 1 1 1 ^ 

3727 - 


' . . . 1 . < . . 1 . 





7300 7350 7400 

Observed Wavelength [A] 

SDSS J0821+3733 

Restframe Wavelength [A] 
3700 3725 3750 



8650 8700 8750 

Observed Wavelength [A] 

SDSS J0830+5116 

Restframe Wavelength [A] 
3700 3725 3750 



8100 8150 8200 

Observed Wavelength [A] 

SDSS J0837+4937 

Restframe Wavelength [A] 
3700 3725 3750 





[O II] 3727 



9000 9050 9100 9150 

Observed Wavelength [A] 

SDSS J0840+5051 

Restframe Wavelength [A] 
3720 3740 4900 5000 



8450 8500 8550 8600 

Observed Wavelength [A] 

SDSS J0841+5017 

Restframe Wavelength [A] 
3700 3725 3750 



8000 8050 8100 8150 

Observed Wavelength [A] 

SDSS J0915-0055 

Restframe Wavelength [A] 
3720 3740 4950 5000 



7 3 
S 2 



L [O II] 3727 J L 



[Oiii 



[Oiii] 4959 




■ [0 11] 


3727 - 


' 1 1 1 1 1 1 1 1 






7150 9300 9400 9500 9600 
Observed Wavelength [A] 

SDSS J0941-0104 



8050 8100 8150 8200 

Observed Wavelength [A] 

SDSS J0944-0147 



7450 7500 9900 100( 

Observed Wavelength [A] 

SDSS J1016-0208 



Fig. 1. — Spectroscopically discovered background emission-lines in BELLS target galaxies. The black solid-line shows the 
BOSS observed flux density, fx, as a function of observed (lower axis) and rest-frame (upper axis) wavelengths. The blue dashed-line shows 
the BOSS template fltted to the continuum of the foreground-galaxy, and the red vertical line show the wavelength of the background 
emission-lines discovered in the differences between the observed flux and the best-flt template. This flgure is continued. 
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Fig. 1. — Continued. Spectroscopically discovered background emission-lines in BELLS target galaxies. 
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Fig. 1. — Continued. Spectroscopically discovered background emission-lines in BELLS target galaxies. 
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Fig. 2. — BELLS Strong Gravitational Lens Galaxies and Foreground Subtracted Images. The 36 Grade-A, B and C lenses 
discovered under HST Cycle 18 program 12209. For each system, the left panel shows the HST ACS-WFC F814W rectified and combined 
images with north up and east to the left, and the right panel reveals the lensed features which remain after the foreground-galaxy has 
been subtracted by the B-spline method described in §3. Probable non-lensed extraneous features that were masked are circled. System 
properties from the BOSS data are listed in Table 2, and those from the HST-ACS data are listed in Table 3. Comments justifying our 
lens grade are provided in Table 4. This figure is continued. 



Einstein radii for higher-redshift lens galaxies. 

These steps were carried out for all BOSS survey plates 
available prior to 2010 June 21: 239 unique plates, with 
an average of approximately 500 successful galaxy red- 
shifts per plate. The systems are listed in Table 2, includ- 
ing the BOSS pipeline redshift, the redshift of our can- 
didate background emission-line source, and the BOSS 
photometry and velocity dispersion measurements. The 
background emission-line spectra remaining after sub- 
tracting the best-fit foreground spectra are presented in 
Figure 1. The original detection of SDSS J1221-0220 
was an artifact of the spectroscopic data-reduction ver- 
sion used at the time of the HST target selection. 

3. IMAGE PROCESSING 

The 45 systems identified in Table 2 were submitted for 
observation under HST Cycle 18 Program ID GO-12209, 
with one single-orbit visit per target using the /-band 
F814W filter, with four sub-exposures of approximately 



500 s each in the ACS-WFC-DITHER-BOX pattern. 

After each complete (and successful) visit, the indi- 
vidual fiat-fielded (FLT) sub-exposure files generated by 
the Space Telescope Science Institute (STScI) CALACS 
pipeline were downloaded, and reduced using ACSPROC 
— a custom-built pipeline that was generalized from the 
SLACS version of the software described in Bolton et al. 
(2008). For each visit, the steps in the current version 
of ACSPROC, implemented in the IDL programming lan- 
guage, are as follows. 

1. Extract a 1500 x 1500 pixel section of the FLT 
file for each sub-exposure, centered on the tar- 
get, performing cosmic-ray rejection using L.A. 
COSMIC (van Dokkum 2001), and sky- level subtrac- 
tion using the FLT MDRIZSKY FITS header, as de- 
termined by the STScI pipeline MULTIDRIZZLE re- 
duction. 

2. Obtain the approximate pixel-shift across sub- 
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Fig. 2. — Continued. BELLS Strong Gravitational Lens Galaxies 
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Fig. 2. — Continued. BELLS Strong Gravitational Lens Galaxies and Foreground Subtracted Images. 
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Fig. 2. — Continued. BELLS Strong Gravitational Lens Galaxies and Foreground Subtracted Images. 
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exposures through fast Fourier transform cross cor- 
relation and compute the centroid of the target 
galaxy for each sub-exposure by fitting an ellipti- 
cal Moffat profile using the MPFIT2DPEAK nonlinear 
least-squares fitting function. 

3. Rectify the sub-exposures onto a uniform 0^/05 per 
pixel grid, with north up and east left, centered 
on the computed centroid. Rectify, with identical 
sampling, an appropriate model point-spread func- 
tion (PSF) which was computed using the Tiny 
Tim EST PSF simulator (Krist 1993). 

4. Combine the sub-exposures into a single stacked 
image, with additional cosmic ray rejection. Simi- 
larly, combine the PSF model. 

5. Within a 280x280 (14''xl4'') section of the stacked 
image centered on the centroid of the target, sepa- 
rately mask potential lensed features, and probable 
non-lensed stars and galaxies. 

6. Fit a set of elliptical radial B-spline models (Bolton 
et al. 2006, 2008) to the masked stacked image us- 
ing the MPFIT nonlinear least-squares implementa- 
tion of the Levenberg-Marquardt algorithm (More 
& Wright 1993), allowing for either none, or com- 
binations of higher-order multipole moments, in- 
cluding the dipole, quadrupole, and octopole mo- 
ments. Each of the B-spline models are manually 
inspected, and the masks are iteratively refined by 
the appearance of strong-lensing features in the B- 
spline-subtracted residuals. The B-spline fit which 
leads to the lowest without introducing obvi- 
ously spurious small-scale structure is selected as 
the model for lens galaxy subtraction. 

Step 4 leads to a combined and rectified image for all of 
the systems listed in Table 2. Step 5 leads to a Grade- A 
foreground-subtracted image for systems with clear and 
convincing evidence of multiple imaging for a subset of 
those systems. Systems with lensed features that are too 
low in surface brightness or contaminated by significant 
bright star light are classified as Grade-B systems and 
systems that are unsuitable for lens modeling due to the 
absence of clearly identifiable counter-images are classi- 
fied as Grade-C systems. The systems without any iden- 
tifiable lensed features are classified as Grade-X systems 
for which we do not complete Step 5. 

The results of the ACSPROC image processing, includ- 
ing the combined and rectified stacked images and the 
foreground-subtracted images are shown in Figure 2 for 
all of the Grade- A, B and C systems. The combined and 
rectified stacked images for the Grade-X systems are in- 
cluded in the Appendix Figure Al for completeness. 

4. PHOTOMETRIC ANALYSIS 

To compute standardized model magnitudes, effective 
radii, projected axis ratios and position angles for the 
BELLS targets, we fit the images to two-dimensional el- 
lipsoidal de Vaucouleurs luminosity profiles. These fits 
are performed over a 51'^ x 51^' square region centered on 
the target galaxies (approximately half the narrower di- 
mension of the EST ACS-WFC CCD aperture in which 
the targets were roughly centered) . The 14^' x 14^' masks 
created with our ACSPROC software described in Step 5 
of §3 are applied in the central regions; stars and neigh- 
boring galaxies outside the manually masked area are 
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Fig. 3. — Distribution of confirmed SLACS and BELLS lens 
HST /814-band magnitudes, shown as a function of lens redshift. 
SLACS lenses are shown with red squares, and BELLS lenses are 
shown with blue diamonds. Stellar masses are indicated by the 
1{)^^Mq upper track and IQ^'^Mq lower track with black dashed- 
lines, based upon the stellar population model of Maraston et al. 
(2009). 




Fig. 4. — Distribution of the ratio between the Einstein radii 
and the effective radii, as a function of lens redshift. SLACS lenses 
are shown with red squares, and BELLS lenses are shown with blue 
diamonds. 



masked with a single-step "clipping" of pixels that de- 
viate by more than four standard deviations from the 
model. The fits are performed using the MPFIT2DFUN 
procedure in the IDL programming language, using a 
flat- weighted inverse variance, and include convolution 
with the appropriate rectified and stacked Tiny Tim 
PSF (Krist 1993). The initial optimization is done by 
sampling the model at one point per data pixel, and re- 
peated with a 5 X 5 sub-sampling per pixel. 

The /814-band apparent magnitudes are computed 
from the full (not truncated) analytic integral of the best- 
fit de Vaucouleurs model, and provided in Table 3. We 
also report the corrections for Galactic extinction from 
Schlegel et al. (1998). Photometric stellar masses are 
computed from the extinction corrected /8i4-band ap- 
parent magnitudes using the Maraston et al. (2009) LRG 
photometric stellar population model with a Salpeter ini- 
tial mass function. We use the intermediate axis effective 
radius (the geometric mean of the major and minor axes). 

Figure 3 shows that the distribution of grade-A lenses 
in HST /814-band magnitudes is similar to that of the 
SLACS sample, lying at higher redshifts, but within the 
same stellar mass range of lO^^M© < < lO^^M©. 
The distribution in the ratio OE/Ref^ of the Einstein ra- 
dius to the lens effective radius as a function of the lens 
redshift is shown in Figure 4. The BELLS sample ex- 
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TABLE 3 

HST-ACS /814-BAND Properties of BELLS Systems 





(1') 




-fi-eff V ) 




PA 

yo) 


1:1 neeifir'n + ir^n 

ICLoolliL/CLtHJll 

\ ' ) 


SDSS J0151+0049 


19.816 


0.049 


0.665 ± 0.002 


0.604 ±0.001 


109.0 ±0.1 


E-S-A 


SDSS J0747+5055 


18.923 


0.111 


1.089 ± 0.002 


0.737 ±0.001 


6.9 ±0.1 


E-S-A 


SDSS J0747+4448 


19.417 


0.066 


0.924 ± 0.002 


0.645 ±0.001 


140.2 ±0.1 


E-S-A 


SDSS J0757+4313 


18.500 


0.068 


3.818 ±0.009 


0.572 ±0.001 


49.8 ±0.1 


L-M-X 


SDSS J0801+4727 


19.911 


0.103 


0.499 ± 0.001 


0.951 ±0.002 


70.8 ± 1.3 


E-S-A 


SDSS J0821+3733 


19.219 


0.079 


0.551 ±0.001 


0.743 ± 0.001 


172.8 ±0.1 


E-S-C 


SDSS J0830+5116 


19.332 


0.087 


0.969 ± 0.002 


0.699 ±0.001 


121.6 ±0.1 


E-S-A 


SDSS J0837+4937 


19.462 


0.060 


0.669 ±0.001 


0.473 ±0.001 


51.3 ±0.1 


L-S-X 


SDSS J0840+5051 


19.984 


0.043 


0.357 ±0.001 


0.699 ±0.001 


63.7 ±0.1 


E-S-C 


SDSS J0841+5017 


19.624 


0.045 


0.648 ± 0.001 


0.939 ± 0.002 


79.3 ±0.8 


E-S-C 


SDSS J0915-0055 


18.518 


0.070 


1.219 ±0.002 


0.925 ±0.001 


173.2 ±0.5 


E-S-B 


SDSS J0941-0104 


20.005 


0.063 


0.458 ± 0.001 


0.568 ±0.001 


129.9 ±0.1 


L-S-C 


SDSS J0944-0147 


19.965 


0.067 


0.478 ± 0.001 


0.785 ±0.002 


91.7 ±0.2 


E-S-A 


SDSS J1016-0208 


19.797 


0.072 


0.465 ± 0.001 


0.816 ±0.002 


1.5 ±0.3 


E-S-C 


SDSS J1039-0014 


19.345 


0.094 


0.812 ±0.001 


0.440 ±0.001 


97.9 ±0.0 


L-S-C 


SDSS J1117-0133 


18.896 


0.079 


2.397 ±0.006 


0.938 ±0.002 


135.3 ± 1.0 


U-U-X 


SDSS J1159-0007 


19.463 


0.049 


0.958 ± 0.002 


0.966 ±0.002 


23.5 ± 1.9 


E-S-A 


SDSS J1215+0047 


19.544 


0.050 


0.651 ±0.001 


0.684 ±0.001 


152.9 ±0.1 


E-S-A 


SDSS J1221-0220 


18.942 


0.055 


0.710 ±0.001 


0.586 ±0.001 


54.9 ±0.1 


E-S-X 


SDSS J1221+3806 


19.984 


0.029 


0.470 ± 0.001 


0.838 ±0.002 


65.6 ±0.3 


E-S-A 


SDSS J1234-0241 


19.269 


0.064 


1.054 ±0.002 


0.762 ±0.002 


93.7 ±0.2 


E-S-A 


SDSS J1318-0104 


19.873 


0.050 


0.687 ±0.002 


0.761 ±0.002 


106.8 ±0.3 


E-S-A 


SDSS J1337+3620 


18.603 


0.023 


2.034 ± 0.003 


0.960 ±0.001 


111.4 ± 1.0 


E-S-A 


SDSS J1344+3258 


19.581 


0.031 


0.524 ±0.001 


0.746 ±0.001 


117.8 ±0.1 


E-S-C 


SDSS J1345-0129 


21.877 


0.071 


1.000 ±0.003 


0.000 ±0.001 


0.0 ±0.1 


E-S-X 


SDSS J1349+3612 


18.555 


0.025 


1.886 ±0.003 


0.743 ±0.001 


8.7 ±0.1 


E-S-A 


SDSS J1352+3216 


19.514 


0.024 


0.579 ± 0.001 


0.949 ±0.001 


127.2 ±0.8 


E-S-A 


SDSS J1452+3323 


19.487 


0.026 


0.623 ± 0.001 


0.836 ±0.001 


87.9 ±0.2 


E-S-B 


SDSS J1503+3225 


20.118 


0.032 


0.769 ± 0.003 


0.625 ±0.002 


82.4 ±0.2 


E-M-B 


SDSS J1522+2910 


19.534 


0.046 


0.890 ± 0.002 


0.579 ±0.001 


130.3 ±0.1 


E-M-A 


SDSS J1537+0220 


19.682 


0.111 


0.386 ± 0.001 


0.694 ±0.001 


52.9 ±0.1 


E-S-C 


SDSS J1541+1812 


19.648 


0.064 


0.759 ± 0.002 


0.755 ±0.002 


151.3 ±0.2 


L-S-A 


SDSS J1542+1629 


18.580 


0.054 


0.726 ± 0.001 


0.786 ±0.001 


89.5 ±0.1 


E-S-A 


SDSS J1545+2748 


18.594 


0.050 


2.589 ± 0.005 


0.661 ±0.001 


105.2 ±0.1 


E-S-A 


SDSS J1601+2138 


19.712 


0.127 


0.436 ± 0.001 


0.960 ±0.001 


110.5 ± 1.0 


E-S-A 


SDSS J1611+1705 


19.611 


0.085 


0.998 ± 0.002 


0.926 ±0.002 


28.4 ±0.9 


L-S-A 


SDSS J1615+2056 


18.193 


0.146 


5.399 ±0.016 


0.899 ±0.002 


145.8 ±0.8 


E-S-X 


SDSS J1618+1930 


19.944 


0.099 


0.429 ± 0.001 


0.513 ±0.001 


9.2 ±0.1 


L-S-X 


SDSS J1627+1910 


19.477 


0.089 


0.360 ± 0.000 


0.661 ±0.001 


101.6 ±0.1 


E-S-X 


SDSS J1631+1854 


17.859 


0.079 


1.434 ±0.001 


0.929 ±0.001 


166.7 ±0.3 


E-S-A 


SDSS J1637+1439 


19.147 


0.093 


1.037 ±0.002 


0.674 ±0.001 


122.9 ±0.1 


L-S-A 


SDSS J2122+0409 


19.501 


0.147 


0.903 ± 0.002 


0.819 ±0.002 


136.6 ±0.3 


E-S-A 


SDSS J2125+0411 


18.271 


0.150 


0.901 ±0.001 


0.689 ±0.001 


111.0 ±0.1 


E-S-A 


SDSS J2303+0037 


18.721 


0.081 


1.349 ± 0.002 


0.760 ±0.001 


84.4 ±0.1 


E-S-A 



Note. — Column 1 provides the SDSS System Name in terms of truncated J2000 R.A. and decl. in the format HHMM±DDMM. Column 2 
provides HST-ACS /8i4-band apparent magnitudes based on fiat-weighted de Vaucouleurs models, and are quoted in the AB system (Oke & Gunn 
1983) without correction for Galactic extinction. Column 3 provides Schlegel et al. (1998) Galactic extinction values that should be subtracted 
from the apparent magnitudes of Column 2 to give dust-corrected magnitudes. Column 4 provides de Vaucouleurs model intermediate axis effective 
radii. Column 5 provides axis ratios of minor to major axes, B/A, for the de Vaucouleurs models. Column 6 provides the major axis position angles 
for the de Vaucouleurs models, measured east from north. Column 7 provides our classification denoting (1) foreground-galaxy morphology, (2) 
foreground-galaxy multiplicity, and (3) status of system as a lens based on available data. Morphology is coded by "E" for early-type (elliptical 
and SO), "L" for late- type (Sa and later), and "U" for unclassified (galaxies that cannot be unambiguously classed as early- or late-type based on 
the HST-ACS /8i4-band data). Multiplicity is coded by "S" for single and "M" for multiple. Lens grades, as justified in Table 4, are coded by "A" 
for systems with clear and convincing evidence of multiple imaging, "B" for systems with strong evidence of multiple imaging but insufficient S/R 
for definite conclusion and/or modeling, "C" for systems with evidence of imaging such as single arcs or likely background features, but without 
detected counter-images, and "X" for all other systems (non-lenses). 

tends the SLAGS sample to a wider range of Oe/RqE 
while still providing substantial overlap with the range 
of ratios seen in the SLAGS sample. 

Each of the HST-ACS /814-band images was visually 
inspected to classify all targets for multiplicity and mor- 
phology. Systems with two or more foreground-galaxies 
of comparable luminosity are classified as "multiple", 
while systems with only a single dominant foreground- 
galaxy are classified as "single". Morphological classi- 
fication is limited to the categories of "early-type" (el- 
liptical and SO), "late- type" (Sa and later spirals), and 
"unclassified" (generally ambiguous between SO and Sa). 
The classification of observed candidates into lenses and 
non-lenses is made by visual examination of the direct 



and B-spline model-subtracted residuals, based on the 
appearance of arcs, rings, and multiple images centered 
on the position of the foreground-galaxy. Our classifica- 
tion is provided in the Table 3, with our justifications for 
BELLS confirmed lenses provided in Table 4. 

5. STRONG LENSING ANALYSIS 

Figure 2 shows the foreground-subtracted images of 
each BELLS target that completed the ACSPROC image 
reduction process described in §3. These are the lensed 
features which remain after the foreground-galaxy has 
been subtracted. This image is passed forward to the lens 
modeling stage as the zeroth extension of a BIZZLE FITS 
file. In addition, the inverse variance of the foreground- 
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TABLE 4 

BELLS Confirmed Gravitational Lenses 



System Name 


Grade 


Justification 


(1) 


(2) 


(3) 



SDSS 


J0151+0049 


A 


SDSS 


J0747+5055 


A 


SDSS 


J0747+4448 


A 


SDSS 


J0801+4727 


A 


SDSS 


J0830+5116 


A 


SDSS 


J0944-0147 


A 


SDSS 


J1159-0007 


A 


SDSS 


J1215+0047 


A 


SDSS 


J1221+3806 


A 


SDSS 


J1234— 0241 


A 


SDSS 


J1318— U1L)4 


A 


SDSS 


J1337+3620 


A 


SDSS 


J1349+3612 


A 


SDSS 


J1352+3216 


A 


SDSS 


J1522+2910 


A 


SDSS 


J1541+1812 


A 


SDSS 


J1542+1629 


A 


SDSS 


J1545+2748 


A 


SDSS 


J1601+2138 


A 


SDSS 


J1611+1705 


A 


SDSS 


J1631+1854 


A 


SDSS 


J1637+1439 


A 


SDSS 


J2122+0409 


A 


SDSS 


J2125+0411 


A 


SDSS 


J2303+0037 


A 



Simple single component source, lensed into a well modeled arc in the southest and a counter-image. 

Image and counter-image are connected in partial Einstein ring, with an additional blob to the west. 

Well modeled 3 + 1 quad with a faint but a definite counter-image, and a compact single source component. 

Gomplex source structure with multiple arcs and counter-arcs. 

Arc with a compact core to the northwest, and with a clear counter-image. 

Well modeled arc and counter-arc, with a matching tail feature. 

Well modeled arc and counter-image with detailed small scale features. 

Well modeled double with detailed feature correspondence between the image and the counter-image. 
Well modeled arc and counter-arc. The features are diffuse but significant. 

Glear image and counter-image. Possible foreground subtraction artifacts affect the counter-image in the 
north due to proximity to the lens center. Fixing the mass axis ratio to the light axis ratio increases model 
Oe hy 0'/03. 

Well modeled compact arc and counter-image. 

Well modeled extended arc with a compact counter-image. 

Well modeled arc and a compact counter-image. Lens galaxy has a slightly disturbed morphology in its core. 
Double image lens with detailed feature correspondence reproduced by the pixelized source plane model. 
Double image lens with feature correspondence. 

Well modeled compact double. Diffuse component to the west reproduced by the pixelized model. 
Well modeled arc and counter-arc. 

Gomplex quad/double transition. The detailed features are reproduced by the pixelized model. 
Well modeled compact image and counter-image. 
Double image lens with feature correspondence. 

Spectacular high surface brightness Einstein ring with multiple source components. 
Well modeled arc and counter-arc. 

Well modeled arc and counter-arc with individual knots that are reproduced by the pixelized source model. 
Extended arc and counter-arc with detailed substructure reproduced by the pixelized source model. 
Glear 3 + 1 quad morphology. There is a galaxy group 10^^ to the east that may contribute shear at the lens 
position and account for the model short comings. 



SDSS J0915— 0055 B Plausibly modeled image and counter-image, including substructure reproduced by the pixelized source 

model. Glassified Grade-B because of the low surface brightness and the significant bright star light con- 
tamination in the south. 

SDSS J1452+3323 B Very faint possible Einstein-ring features, but the surface brightness and S/R are too low for modeling. 
SDSS J1503+3225 B Very plausible image-counter-image pair to the northwest and southeast, with detailed feature correspon- 
dence, but the group environment of the lens complicates modeling and interpretation. 



SDSS J0821+3733 G Likely background features seen in the residual image, but no counter-images. 

SDSS J0840+5051 G Likely background features seen in the residual image, but no counter-images. 

SDSS J0841+5017 G Extended source to south inconsistent with the compact source to the north as image-counter-image pair. 

SDSS J0941— 0104 G Glear distorted arc, but no detected counter-image. 

SDSS J1016— 0208 G Tangentially elongated image to the northeast with no detected counter-image. 

SDSS J1039— 0014 G Disky late-type structure with lensed arc to the east, but no apparent counter-image. 

SDSS J1344+3258 G Glear evidence of background-galaxy to the west, but no evidence of a counter-image. 

SDSS J1537+0220 G Pronounced asymmetry in the lens galaxy brightness profile; plausible background-galaxy arc to the south- 
west, but no evidence of a counter-image. 



Note. — Column 1 provides the SDSS System Name in terms of truncated J2000 R.A. and decl. in the format HHMM+DDMM. Column 2 
provides our lens-grade, coded by "A" for the 25 systems with clear and convincing evidence of multiple imaging, "B" for the 3 systems with strong 
evidence of multiple imaging but insufficient S/R for definite conclusion and/or modeling, "C" for the 8 systems with evidence of imaging such as 
single arcs or likely background features, but without detected counter-images. Column 3 provides our lens-grade justification. SIE model parameters 
for the grade- A lenses are listed in Table 5. The remainder of the observed sample, listed in Table Al, receives a grade of "X" (non-lenses). 



subtracted image, the lensed features mask, the mask for 
extraneous unlensed features, and the HST ACS-WFC 
PSF are included in successive FITS extensions. The 
best-fit B-spline model parameters of the foreground- 
galaxy are included in the FITS header so that the cen- 
ter of the lens model may be constrained by the lumi- 
nous distribution, and so that the ellipticity of the vis- 
ible galaxy may be used to initialize the lens model. 
These BELLS BIZZLE files are input into our ACSLENS 
software, implemented in the Python programming lan- 
guage, which accommodates the SIE lens model (Kas- 
siola & Kovner 1993; Kormann et al. 1994; Keeton & 
Kochanek 1998), and both parametric or pixelized source 
models. 

The SIE model is a two-dimensional potential of simi- 
lar concentric and aligned elliptical isodensity contours, 
with axis ratio qsiE- The projected surface density of the 



SIE falls off as 1/i?, and the model provides flat galaxy 
rotation curves. The model is parameterized by its angu- 
lar Einstein radius, Oe^ which is related to the physical 
mass model through 



47r- 



'SIE 



Dls 



Ds 



(1) 



where asiE is a velocity dispersion parameter and Dls 
and Ds are cosmological angular diameter distances from 
lens and observer to source, respectively. We adopt 
the intermediate- axis normalization of Kormann et al. 
(1994), whereby the mass within a given isodensity con- 
tour remains constant at fixed Oe for changing gsiE- 

We attempt to model all candidate lenses with an SIE 
lens-mass model and a parametric source-plane surface- 
brightness model, where the source may be composed of 
either single or multiple Gaussian, de Vaucouleurs (1948) 
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Fig. 5. — SIE Lens Models with Sersic Source Models as compared to Pixelized Source Models. SIE Lens Models are shown 
for each of the 25 confirmed-to-date Grade-A lenses discovered under HST Cycle 18 Program ID GO-12209. For each system, the left panel 
shows the lensed features which remain after foreground-galaxy subtraction by the B-spline method described in §3, repeated from the 
rightmost panel of Figure 2. Best-fit SIE Models with Sersic sources are shown in the second panel from the left, including the residuals 
of the data minus the parametric model in the third panel, with the de-lensed Sersic source model shown in the fourth panel, at a scale 
optimized for the source plane. Best-fit SIE Model parameters are provided in Table 5. SIE Models with pixelized sources are shown in 
the fifth panel from the left, including the residuals of the data minus the pizelized model in the sixth panel, with the de-lensed pixelized 
source model shown in the rightmost panel, at the same scale as used for the parametric source model. Comments justifying our lens grade 
are provided in Table 4. The color scale is provided in Figure 2. This figure is continued. 
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Fig. 5. — Continued. SIE Lens Models with Sersic Source Models as compared to Pixelized Source Models. 
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TABLE 5 

BELLS Grade-A Strong Lens SIE Model Parameters 



System Name 




Q'SIE 


P.A. (°) 


Ns 


m8i4 


11 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


SDSS J0151+0049 


0.676 


0.752 


111.0 


1 


22.51 


8.71 


SDSS J0747+5055 


0.754 


0.641 


4.9 


2 


21.46 


2.95 


SDSS J0747+4448 


0.610 


0.723 


147.1 


1 


23.77 


39.72 


SDSS J0801+4727 


0.491 


0.891 


41.1 


2 


22.07 


3.82 


SDSS J0830+5116 


1.142 


0.887 


137.0 


2 


21.89 


7.43 


SDSS J0944-0147 


0.725 


0.922 


108.4 


1 


23.29 


5.09 


SDSS J1159-0007 


0.683 


0.815 


125.3 


2 


21.92 


9.46 


SDSS J1215+0047 


1.368 


0.742 


123.2 


1 


21.25 


3.69 


SDSS J1221+3806 


0.699 


0.745 


91.4 


3 


22.50 


3.43 


SDSS J1234-0241 


0.533 


0.279 


77.7 


2 


21.92 


2.56 


SDSS J1318-0104 


0.679 


0.836 


105.2 


1 


22.87 


6.13 


SDSS J1337+3620 


1.386 


0.682 


137.4 


2 


21.71 


11.55 


SDSS J1349+3612 


0.750 


0.711 


177.5 


1 


22.06 


4.79 


SDSS J1352+3216 


1.823 


0.859 


115.2 


1 


21.97 


5.65 


SDSS J1522+2910 


0.736 


0.741 


143.7 


1 


22.24 


5.55 


SDSS J1541+1812 


0.640 


0.927 


142.4 


1 


23.80 


9.76 


SDSS J1542+1629 


1.042 


0.806 


105.6 


1 


22.32 


3.28 


SDSS J1545+2748 


1.210 


0.418 


99.0 


2 


22.58 


4.35 


SDSS J1601+2138 


0.858 


0.906 


164.6 


1 


22.84 


3.50 


SDSS J1611+1705 


0.580 


0.740 


4.3 


2 


22.13 


2.27 


SDSS J1631+1854 


1.634 


0.878 


35.3 


3 


20.41 


19.15 


SDSS J1637+1439 


0.650 


0.750 


122.5 


1 


22.99 


10.65 


SDSS J2122+0409 


1.580 


0.634 


119.7 


2 


22.31 


5.73 


SDSS J2125+0411 


1.204 


0.821 


83.0 


1 


21.11 


4.32 


SDSS J2303+0037 


1.016 


0.391 


88.8 


1 


23.10 


8.06 



Note. — The method used to compute the best-fit SIE model pa- 
rameters is described in §5. Column 1 provides the SDSS System 
Name in terms of truncated J2000 R.A. and decl. in the format 
HHMMiDDMM. Column 2 provides the best-fit SIE model Einstein 
radius. Column 3 provides the best-fit minor to major axis ratio for the 
SIE models. Column 4 provides the major axis position angle for the 
SIE models, measured east from north. Column 5 provides the number 
of Sersic models added to the source plane to provide the best-fit. Col- 
umn 6 provides the /8i4-band magnitude of the lensed source. Column 
7 provides the magnification of the lensed source. 

or Sersic (1968) distributions, as necessary to obtain a 
good (x^) fit. The model lensed image is generated by 
ray-tracing through the analytic SIE mass model, and 
then convolved with the EST ACS-WFC PSF. After 
choosing reasonable initial parameter values, the SIE and 
source-plane model parameters are optimized nonlinearly 
with the least sq Levenberg-Marquardt implementation 
in the scipy. optimize package. In the SIE mass-model 
optimization, the mass centroid is fixed to the centroid 
of the luminous distribution as determined from the B- 
spline lens-galaxy photometric analysis, while the SIE 
strength, axis ratio, and major-axis position angle are 
all allowed to vary. The final outcome is a set of lens- 
model and source-component parameters, along with a 
model for the lensed images. 

Our best-fit parametric source model is then comple- 
mented with a pixelized source-plane model, which we 
compute through a matrix inversion including linear reg- 
ularization, without re-optimizing the lens model (War- 
ren & Dye 2003; Treu & Koopmans 2004; Wayth & Web- 
ster 2006; Koopmans et al. 2006). 

If necessary to obtain a better correspondence between 
the parametric and pixelized model images, the appar- 
ent multiplicity and spatial orientation of the pixelized 
source distribution are used to place additional Sersic 
sources, and then the parametric model is re-optimized 
together with the SIE mass-model parameters. Several 
iterations are pursued until the parametric and pixelized 
models are in reasonable agreement. Although we do not 
use the pixelized models in the lens model optimization, 
their ability to capture more detail in the lensed galax- 



ies serves as an important confirmation of the lensing 
hypothesis in most cases. 

Our best-fit SIE lens models with Sersic source models, 
along with our pixelized source models, are presented in 
Figure 5, and the best-fit model parameters are listed 
in Table 5. Our justifications for the BELLS confirmed 
lenses grading are provided in Table 4. 

6. SUMMARY AND CONCLUSIONS 

We present 44 of the 45 targets observed-to-date in 
the BOSS Emission-Line Lens Survey (BELLS) from the 
EST ACS-WFC Cycle 18 Program ID GO-12209 imag- 
ing program. The BELLS catalog presently includes 25 
grade-A strong galaxy-galaxy lenses complete with lens 
and source redshifts, /814-band foreground-galaxy pho- 
tometry, and SIE models. This represents a significant 
expansion of the SLACS survey, extending the lens red- 
shifts to — 0.66, and with source redshifts extended 
to — 1.52. The combined SLACS and BELLS catalog 
provides a homogeneous sample of strong lens galaxies at 
significant cosmological look-back times. Future papers 
in this series will apply the combined SLACS and BELLS 
samples to the analysis of the structure and evolution of 
massive elliptical galaxies. 

We have described the use of our ACSPROC RST ACS- 
WFC image processing software to make luminosity mea- 
surements and to produce our BIZZLE data product. We 
then analyze the BIZZLE data using our ACSLENS lens- 
modeling software to make parameterized and pixelized 
lens models. Our analysis demonstrates that simple 
SIE lens models, combined with multiple Sersic ellipsoid 
models of the lensed background-galaxies, can reproduce 
the lensed image configurations in great detail, and our 
pixelized source models reproduce more detailed features 
in the lensed images and add finer structure to the de- 
lensed source model. 

The high success rate of our HST imaging program 
demonstrates that the spectroscopic strong-lens selection 
methodology can be successfully applied to the BOSS 
spectroscopic database even though BOSS is targeting 
a more distant — and hence fainter — galaxy popula- 
tion than SDSS-I. The probability of finding a galaxy- 
galaxy lens using the spectroscopic discovery method was 
considered in Dobler et al. (2008) for the SLACS sur- 
vey, including selection effects such as the finite size of 
the spectroscopic fiber (which selects against large sep- 
aration lenses), and the effectiveness of spectral noise 
modeling (which selects against sources that have red- 
shifted emission-lines coincident with strong emission- 
lines in the sky). The probability that a galaxy-galaxy 
lens would be spectroscopically discovered from either 
the SDSS-I or SDSS-III spectra, and lead to a confirmed 
lens, is under consideration in Arneson et al. (2011) using 
a Monte Carlo simulation to measure the effects of varia- 
tion in the lens mass, ellipticity and morphology, source 
extent, and the logarithm slope of the lens luminosity 
profile. 

The strong lensing measurements presented in this 
work afford a unique opportunity to test the results of 
theoretical models and numerical simulations of galaxy 
formation, merging, and evolution. This is due to the fact 
that strong lensing measures the projected mass inside 
the aperture defined by the Einstein radius, in a nearly 
model-independent sense, and without need for modeling 
of stellar populations and luminosity evolution. 

In comparison with the SDSS-I spectra of the SLACS 
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lens sample, the BOSS spectra for the BELLS sample 
have significantly lower signal to noise. Hence, the re- 
sulting velocity dispersion measurements for the BELLS 
sample are correspondingly noisier and more care must 
be taken in interpreting these measurements in the con- 
text of dynamical galaxy modeling. Work in prepara- 
tion will present a joint analysis of all BELLS lenses 
that marginalizes over the full velocity dispersion likeli- 
hood function for individual spectra (see Shu et al. 2011). 
Deeper spectroscopy with large ground-based telescopes 
can provide more precise stellar velocity dispersion mea- 
surements. 

We note that the BOSS spectroscopic database as of 
summer 2011, which will constitute the ninth SDSS pub- 
lic Data Release (DR9), contains 819 unique BOSS plates 
with an average of approximately 600 successfully red- 
shifted galaxies per plate. (The increase in the num- 
ber of spectroscopically confirmed galaxies per plate was 
produced by a combination of improvements in target- 
ing, hardware, operations, and analysis software during 
the BOSS commissioning phase and early operations.) 
Thus, DR9 is at least a fourfold increase over the sample 
we used here. The final BOSS sample, to be completed 
in mid-2014 and released at the end of that year, will 
be an order of magnitude increase over the sample used 
here. In addition, ongoing BOSS spectroscopic pipeline 
development (including implementation of the two di- 
mensional PSF-based extraction algorithm of Bolton & 
Schlegel 2010) continues to increase the quality of the 
reduced spectra and redshift /classification analysis for 
BOSS significantly. Hence, there will be at least several 
hundred strong lenses to be found in BOSS spectroscopy 
at the end of the survey, offering the possibility of mea- 
suring the cosmic evolution of the relationship between 



mass-density structure and stellar mass, luminosity, and 
rest-frame color, in great quantitative detail. 
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Fig. Al. BELLS Grade-X Non-Lens Galaxies. The 8 non-lenses, discovered to-date under HST Cycle 18 Program ID GO-12209, 
Table 2 presents their properties as derived from the BOSS data, and Table 3 presents those derived from the HST data. For each system, 
the panel shows the HST ACS-WFC F814W rectified and combined images with north up and east to the left. Comments justifying the 
lens- grade of "X" are provided in Table Al. The color scale is provided in Figure 2. 



TABLE Al 
BELLS Grade-X Non-Lens Galaxies 



System Name 


Justification 


(1) 


(2) 



SDSS J0757+4313 
SDSS J0837+4937 
SDSS J1117-0133 

SDSS J1221-0220 
SDSS J1345-0129 
SDSS J1615+2056 
SDSS J1618+1930 

SDSS J1627+1910 



Two components with spiral structure; any lensing features are difficult to disentangle. 
Late-type with faint candidate features, but with no evidence of multiple imaging. 

Pronounced asymmetry in the foreground-galaxy plus possible spiral and dust features. Possible lensed background 
sources cannot be unambiguously modeled given single-band imaging. 

Double nucleus in the foreground-galaxy, with no evidence of background-galaxy features. 

Faint possible background-galaxy features to the west, but without plausible evidence of strong lensing. 

Possible candidate-arc extending to the northwest and to the east, but the system is difficult to interpret as a lens. 

No obvious background features, but spiral structure and/or dust in the foreground-galaxy complicates 

interpretation. 

Possible faint background-galaxy features to the northeast, but without evidence of strong lensing. 



Note. — Column 1 provides the SDSS System Name in terms of truncated J2000 R.A. and decl. in the format HHMMiDDMM. Column 2 
provides our justification for the lens-grade of "X" (non-lenses). 
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